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Experimental implementation of a quantum algorithm requires unitary operator decomposition. 
Here we treat the unitary operator decomposition as an optimization problem and use Genetic 
Algorithm, a global optimization method inspired by nature's evolutionary process for operator 
decomposition. As an application, we apply this to NMR Quantum Information Processing and 
find a probabilistic way of doing universal quantum computation using global hard pulses. We also 
demonstrate efficient creation of singlet state (as a special case of Bell state) directly from thermal 
equilibrium using an optimum sequence of pulses. 



I. INTRODUCTION 

Quantum computation (QC) may possibly be the most 
remarkable proposal of practical application of quantum 
mechanics [T]. Quantum information processor (QIP) 
is one which can control a large quantum system well 
enough to perform an arbitrary quantum algorithm. It 
holds out tremendous promise for efficiently solving some 
of the most difficult problems in computational science, 
such as integer factorization^ , database search [3J and 
quantum simulation that are intractable on any present 
or future conventional computer [4]-[7]. 

Genetic algorithms (GA) are stochastic global search 
method based on the mechanics of natural biological evo- 
lution [8]. It was first proposed by John Holland in 
1975 [9]. GA operates on a population of solutions of 
a specific problem by encoding the solutions to a simple 
chromosome-like data structure, and applies recombina- 
tion operators. At each generation, a new population 
is created by breeding individuals (selected according to 
their fitness value) together using operators borrowed 
from natural genetics. This process leads to the evo- 
lution of individuals and generate populations that are 
better suited to their environment. GAs are attractive 
in engineering design and applications because they are 
easy to use and are likely to find the globally best de- 
sign or solution, which is superior to any other design or 
solution [ID] . 

Long-lasting coherence and high fidelity controls in nu- 
clear magnetic resonance (NMR) spectroscopy are ideal 
for quantum information processing. It made NMR QIP 
to perform first experimental implementations of ba- 
sic quantum computation algorithms like Deutsch-Jozsa 
[IT] algorithm, Grover's algorithm [3] , preparation of the 
three-qubit GHZ state [12] . Algorithms like shor's factor- 
ization algorithm [2J, dense coding[T3], quantum telepor- 
tation [T3] are demonstrated in NMR QIP and is one of 
the best candidate for testing basic principles of quantum 
mechanics [El [16] . 
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Decomposing a unitary operator in terms of experi- 
mentally preferable operators is the main task in exper- 
imental implementation of a quantum algorithm. There 
are several proposals for doing this for NMR QIP, 
such as SMPS (Strongly Modulated Pulses) by Cory et 
<zL[T7], GRAPE (Gradient Ascent Pulse Engineering) by 
Khaneja et aL|18| and algorithmic approach by Ajoy et 
al. [TDJ • Here we investigate the use of Genetic algorithm 
for direct numerical optimization of pulse sequences and 
devise a probabilistic method for doing universal quan- 
tum computation by using hard pulses. We also inves- 
tigate quantum state preparation using GA optimiza- 
tion. We demonstrate singlet state preparation (along 
with preparation of other three Bell states) directly from 
thermal equilibrium in two qubit homonulear NMR sys- 
tem by using global hard pulses and pulsed-field gradients 
gDJ. InSec|H]ofthis paper, we describe the optimization 
procedure and in Sec III we outline the experimental im- 
plementations. 



II. GENETIC ALGORITHM FOR NMR PULSE 
SEQUENCE GENERATION 

In liquid state NMR, the system Hamiltonian is com- 
posed of the interactions of spins with external magnetic 
field (chemical shifts) and coupling interactions among 
spins. Combining this with external RF pulses (with spe- 
cific frequency, amplitude and phase) can simulate any 
preferred effective Hamiltonian |21j . Hence unitary op- 
erator decomposition problem in NMR can be treated as 
an optimization problem which can give optimal values 
of pulse parameters and delay durations. Here optimality 
is determined by a proper fitness function which depends 
on target Hamiltonian or target state. 

We have performed pulse sequence optimization us- 
ing GA for quantum logic gates (operator optimiza- 
tion) and quantum state preparation (state-to-state 
optimization)^?]- State-to-state optimization converges 
faster than operator optimization (there can be many 
operators which can perform same state-to-state trans- 
fer). In the discussion given below, Single Qubit Rota- 
tion (SQR) pulses and CNOT gates are operator opti- 
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mization whereas creation of pseudo pure state (PPS) 
and Bell states are state-to-state optimization. 
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A. Representation Scheme: 

Representation scheme is the method used for encod- 
ing the solution of the problem to individuals in genetic 
evolution. Designing a good genetic representation that 
is expressive and evolvable is a hard problem in evolu- 
tionary computation [8]. Defining proper representation 
scheme is the first step in genetic algorithm optimization. 

In our representation scheme we have selected the gene 
as a combination of (i) an array of pulses which will ap- 
ply simultaneously on each channel with arbitrary ampli- 
tude (9) and phase (0) and (ii) an arbitrary length delay 
(d). It can be easily shown that repeated application of 
the above gene forms the most general pulse sequence in 
NMR. The individual, which represents a valid solution 
can be represented as a matrix of size (n+1) X 2m and is 
shown in Eqnjl] Here m is the number of genes in each 
individual and n is the number of channels or spins. So 
the problem is to find an optimized matrix (Eqn[lJ where 
optimality condition is posed by a fitness function (Sec. 
TlB I. 
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Iii the beginning of optimization, we select the number 
of genes (m) as an arbitrary number which is a function 
of the complexity of the problem to be solved. If the 
fidelity (fitness of the best individu al in the population, 
calculated using fitness function Sec II B ) crosses a cutoff 



value of (say > 99%), the optimization program tries to 
reduce the gene number by assigning zero value to gene 
parameters, and if not the program will rerun with more 
number of genes. In the optimization, we have used a 
population size of 100 individuals for 1000 generations. 
All the programs are written in Matlab® in combination 
with Matlab's optimization toolbox. 



B. Fitness Function: 



FIG. 1. A twospin system with chemical shifts ±<5 and cou- 
pling J 



the optimization has to run for maximizing the fitness 
function. 



F P ui = Trace(U pu i x u} ar ) 



(2) 



In state-to-state optimization, the optimization pro- 
gram will run over different possibilities of U pu i to pre- 
pare a preferred target state (ptar) from an initial state 
(pin). The fitness function we have selected here is given 
by, 



Fstate = Trace(U pu i x p in x {J p J x p\ ar ) (3) 

Here also the optimization has to run for maximizing 
the fitness function. 



III. TWO QUBIT HOMONUCLEAR CASE: 

Consider a two qubit NMR homonuclear system 
(FigjTJ) with chemical shifts ±5 and coupling J. Assum- 
ing weak coupling approximation (6 3> J), the Hamilto- 
nians can be written as!21l. 



H = H cs + Hj = 2it5(l\ - l2 z) + ^J(lUz) ( 4 ) 

For single qubit rotation (SQR) in the above system, 
one can use spin selective pulses (low power, long R.F 
pulses) which will excite a small spectral region around 
the spin to be selected [22]. On the other hand, by using 
the natural chemical shift difference between two spins, 
we show here how to implement SQR with global hard 
(non-selective) pulses. Later we extend this for perform- 
ing two qubit homonuclear universal quantum computa- 
tion using global hard pulses only. 



A fitness function is a particular type of objective func- 
tion that prescribes the optimality of a solution or indi- 
vidual. In operator optimization, GA tries to reach a 
preferred target unitary operator (Utar) from an initial 
random guess pulse sequence operator (U pu i). Here we 
selected a fitness function which is proportional to the 
projection of the U pu i into the target operator (U tar )- It 
will give a maximum value of 1.0 when U pu i=Utar and 



A. Operator Optimisation 

Operator optimization deals with Pulse sequence gen- 
eration for quantum logic gates. Here we show two 
essential unitary operators for universal quantum com- 
putation: Single Qubit Rotation (SQR) and two qubit 
Controlled-NOT gate. 
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FIG. 2. Pulse sequence for single qubit rotation. First two 
filled pulses are (n/2). Flip angle of the third pulse is (0/2) 
with phase <j>. The brackets above each pulse contain the flip 
angle (first number) and the phase (second number). 



Spin to be Excited 


<h 


4>2 


1 ( 


^ - tt/2) ( 


<t> + tt/2) 


2 ( 


+ 7T/2) ( 


0-7T/2) 



TABLE I. Spin selection and Phase i 
by the Phases cj>i and ^2 



> of the SQR is controlled 



Single Qubit Rotations using non-selective pulses: 
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For a two qubit homonuclear NMR system, we first 
consider the case J — (EqnjSJ). 



H = H CS = 2nS(ll - II) 



(5) 



Evolution under such Hamiltonian can create a rela- 
tive phase among spins. The relative phase accumulated 
is proportional to the chemical shift difference (25) and 
evolution time. Combining this relative phase with global 
rotation hard pulses, single qubit operations can be per- 
formed . 

The optimized pulse sequence for SQR is shown in 
Figj2] In the be ginning we had selected m = 3, i.e. three 
hard pulses and three delays. The optimised sequence 
has three hard pulses and single delay. The flip angle 
(9) of the SQR pulse is determined by the delay and the 
flip angle of the third pulse, whereas the phase of SQR 
(4>) and spin selection is determined by phases of all the 
three pulses (Tabjl]). 

Experimental verification of SQR pulse sequence in 5- 
Bromofuroic acid (Fig ]3(a) I (in CqD 6 ) system is shown 



in Fig 3(c) and 3(d) The total length of the pulse se- 
quence for (7r/2) SQR pulse is less than 500/xs where as 
conventional method (using a selective soft pulse) would 
need a 2ms shaped pulse. This shortening in time can 
lead to a significant improvement in a quantum circuits. 



(d) 



FIG. 3. (a). Chemical structure of 5-Bromofuroic acid. 
Diagonal elements in the table contains the chemical shifts of 
protons at 500 MHz and the nondiagonal element represents J 
coupling (sample dissolved in CeDe). (b). Equilibrium Spec- 
trum (c). (71-/2)2, SQR pulse on spin-I (d). (tt/2) v SQR pulse 
on spin-2. The average experimental fidelity (calculated using 
standard definition |17j ; spectral intensities were compared 
with equilibrium spectrum) for the SQR pulse is obtained as 
99.9%. 




FIG. 4. Matlab simulation study of Fidelity variation with 
different values (J/S) and flip angle (0). 



The above analysis also holds for J / 0, so long as 
7-Bi 3> 6, J; except that introduction of J coupling de- 
phases the final state and results in a fidelity loss. The 
fidelity of the pulse sequence (fitness function Eqn[2]) is 
studied using Matlab simulation (FigQ and is > 99.8% 
for J/S < 0.1 and < tt/2. (J/5 < 0.1 is the limit for 



weakly coupled spins and in this paper we are dealing 
with only weakly coupled spins). For 7r/2 < 8 < tt, the 
theoretical fidelity can reach upto 99.5% (FigQ. 
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2. Controlled NOT gate (C-NOT) 

The Controlled-NOT gate is an essential component 
in the construction of a quantum computer. Any quan- 
tum circuit can be simulated to an arbitrary degree of 
precision using a combination of C-NOT gates and sin- 
gle qubit rotations pQ. The optimised pulse sequence for 
Controlled-NOT is shown in Fig 5(a) (obtained using the 
Hamiltonian of Eqn. |4|) 
obtained by tuning 



All the four C-NOT gates can be 
as shown in Tab pi} The pulse 
sequence is identical for all the four C-NOT gates except 
the angles 6 and <fr. 



X3ir/4) (»,*) (7w/4,y) (w/4,5) (0,y) 



D 



0.5 



1 

4 5 



21 



+ 6 



(a) 



mm 

-0.5 ■— 1 -0.S 1 -0.5' — -0.5' ■ ■ ■ ■ 

1234 1234 1234 1234 1234 

CO (") (iH) (vi) (v) 

(b) 




(c) 



FIG. 5. Pulse sequence for Controlled-NOT gate. The brack- 
ets above each pulse contain the flip angle (first number) and 
the phase (second number), (b) Diagonal element tomogra- 
phy of (i). Equilibrium state (ii). after applying C-NOT(l,2) 
(iii). C-NOT(2,l), (iv). C-NOT (1,2) and (v). C-NOT (2,1). 
Here {1, 2, 3, 4} corresponds to {jOO), |01), |10), |11)}. An 
average experimental fidelity of 99.9% is observed, (c) The 
fidelity (F) versus (J/S) plot of C-NOT gate. 

The experimental implementation of various C-NOT 
gates in 5-Bromofuroic acid is shown in Fig |5(b)| An 
average experimental fidelity (calculated using 17] with 
only diagonal element of the density matrix) of 99.9% is 
achieved. 

The theoretical fidelity of the operator (using Eqn{2| is 
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TABLE II. 



values for all four CNOT gates 



however dependent on the ratio (J/S). Matlab simulation 
of fidelity with (J/S) for C-NOT gate pulse sequence is 
shown in Fig |5fcJ| and is > 99.99% for (J/S) = 0.01 and 
> 99.84% for (J/S) < 0.1. This means that even if one 
needs 10 C-NOT gates in a quantum circuit, the fidelity 
is >99%. 



B. State to State Optimization 

State to state optimization deals with pulse sequence 
generation for quantum state preparation. Here we 
added gradient pulses to the optimization procedure, 
which enables to perform non-unitary transformations 
and show two important quantum state preparations: 
Pseudo Pure state creation and Bell state creation di- 
rectly from thermal equilibrium (mixed) state using 
global hard pulses. 



1. PPS creation 

Quantum information processing by NMR spec- 
troscopy uses pseudo-pure states to mimic the evolution 
and observations on pure states [23]. There are several 
methods for creating pseudo pure states (PPS) from ther- 
mal equilibrium [23-25J. Here we closely follow the Spa- 
tial Averaging Method (SAM) of Cory et al. However, 
the SAM uses spin selective pulses, which in homonuclcar 
system becomes soft long pulses. Here we obtain a novel 
pulse sequence using only non-selective (hard) pulses for 
homonuclear two qubit system (Eqnj4|. The optimiza- 
tion problem here is a state to state optimization with 
thermal equilibrium state Ap eq = 11+1% as the initial 
state and Ap 00 = J* + + 2I*J| as the final state. 

For easier optimization (and experimental implemen- 
tation) we have fixed all the pulses to be (7r/2) and op- 
timized only the pulse phases. The sequence consists of 
six (tt/2) pulses and one (tt) pulse for refocusing chem- 
ical shift (Fig 6(a)). The phase of the (tt) pulse can be 
controlled to achieve either 1 00) pps or |11) pps. The 



other PPS are obtained by using a combination sequence 
of PPS and SQR 7r pulse. The experimental results are 
shown Fig 6(b) An average experimental fidelity (calcu- 
lated using [17j with only diagonal element of the density 
matrix) of 99.5% is obtained for various PPS. 

The theoretical fidelity (using Eqn[3| of PPS prepara- 
tion pulse sequence is also dependent dependent on the 
ratio J / S and studied using Matlab simulation (Fig 6(c)[ ) . 
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FIG. 6. (a)Pulse sequence for PPS creation. All filled pulses 
are {tt/2) and nonfilled are ir with phases written above it. ± 
in the phase of n pulse determines the pps to be created ( 1 00) 
or 1 11}). (b) Tomography of diagonal elements after preparing 
(i) |00) pps, (ii) |01) pps, (iii) jlO) and (iv) |11) pps. Here {1, 
2, 3, 4} corresponds to {|00), |01), |10), |11)}. An average 
experimental fidelity of 99.5% was obtained, (c) The fidelity 
(F) — (J/S) plot of PPS generation pulse sequence (Fig 6(a) I. 



The theoritical fidelity of the pulse sequence is > 99.9% 
for J/S < 0.1. 



decided by controlling the phase of the pulses and delay 



durations according to Tab III 
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TABLE III. The values of cj> and d in the pulse sequence shown 



in Fig 7(a) 



The experimental implementation is performed in 5- 
Bromofuroic acid, for the Bell state \<f>~) = ^(|01) — 
1 10)), also known as the singlet state and is a long lived 
state for the Hamiltonian H = J1.J2 [2"9"I I3T)] . 
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The experimental results are shown in Fig |7(b)1 and 



7(c) Since p^- is not directly observable, we convert 



the created singlet state into observable single quantum 



coherence by applying U 



x e 



-i(7T/4)(7i-lJ) 



An experimental fidelity (using the standard definition in 
|17j ) more than of 99.5% is achieved. This is the shortest 
known pulse sequence for creating pure singlet state in a 
two qubit homonuclear NMR system. 

The singlet state life time (T s ) has been measured by 
applying Waltz-16 [3T] spin lock sequence for a variable 
time period (0 — 20 Sec) and is obtained as 11.2 Sec 
(Figj§ longer than T x = 8.7 Sec and T 2 = 3.8 Sec. 

Pulse sequence generation using GA for higher qubit 
systems is in progress. SQR and C-NOT gates using hard 
pulses (Sec. IIIAl and III A 2 1 is valid for higher spin 
systems with homonuclear spin pairs (for example H — 1 
H - 19 F - 19 F system in 2,3-Difluro-6-nitrophenol[2"2"]). 



2. Creation of Bell states directly from equilibrium state 



IV. CONCLUSION 



Bell states are maximally entangled two qubit states 
(also known as Einstein- Podolsky- Rosen states) [25] . 
They play a crucial role in several applications of 
quantum computation and quantum information theory. 
They have been used for teleportation, dense coding and 
entanglement swapping [T31 HU [57]. Creation of Bell 
state is more demanding in quantum computation and 
conventionally (in NMR) it requires product state pps 
creation + Hadamard gate + C-NOT gate [55] ■ Here we 
integrated all these in a single pulse sequence and opti- 
mize with GA. Again we have kept all pulse amplitude 
to be (7r/2) and optimized for pulse phases and delay 
durations. The optimized pulse sequence (Fig 7(a) ) has 
only ten non-selective pulses. The final Bell state can be 



In conclusion, we have used the global optimization 
power of Genetic algorithm for (i) efficiently implement- 
ing SQR and C-NOT gates and (ii) creating PPS in 
homonuclear two qubit system using only hard pulses, 
and hence showed a method for doing universal quan- 
tum computation in such systems. We have also demon- 
strated the creation of Bell states directly from thermal 
equilibrium state. The pulse sequence for Bell state cre- 
ation is the shortest known sequence. Unlike GRAPE 
[IB] or SMPS [T7J, all the pulse sequences discussed here 
are valid for any weakly coupled homonuclear spin pairs. 
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FIG. 7. (a). Pulse sequence for creating Bell states directly form thermal equilibrium state. All filled pulses are (tt/2) and 
nonfilled are n. The phase of each pulse is written above. The values of '0's and 'd's are shown in Tab jIII| (b). Observing singlet 

state after applying U(= e~ i( - 7v/2)I ^ 



x e '("V 4 )^* r *^), (c). Density matrix tomography of created singlet state. 




FIG. 8. Anti-phase signal decay as a function of interval and 
fits to a single exponential decay. The initial amplitude of 
singlet state is normalized to one. Observed Singlet state life 
time is T s — 11.2 Sec. The system has a Ti = 8.7 Sec and 
T 2 = 3.8 Sec. 
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